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Abstract 

We state and prove several theorems that demonstrate how the coordinate Bethe Ansatz 
for the eigenvectors of suitable transfer matrices of a generalised inhomogeneous five- vertex 
model on the square lattice, given certain conditions hold, is equivalent to the Gessel-Viennot 
determinant for the number of configurations of N non-intersecting directed lattice paths, 
or vicious walkers, with various boundary conditions. Our theorems are sufficiently general 
to allow generalisation to any regular planar lattice. 
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1 Introduction 



The problem of non-intersecting paths or vicious walkers has been studied by the statistical 
mechanics community who have been interested in them as simple models of various polymer 
and other physical systems and independently by the combinatorics community who have been 
interested in them in connection with binomial determinants. In statistical mechanics they 
are generally known as vicious walkers, a term coined by Fisher Q, who studied the continuous 
version of the model. Various cases of the lattice problem were later studied by Forrester |5], ^, [?]] . 
Independently in the area of combinatorics the problem of non-intersecting paths was solved by 
a very general theorem of Gessel and Viennot |)|, ||, following the work of Lindstrom [14], and 



Karlin and McGregor [12, 11 1. All these studies express the number of configurations as the 
value of a determinant. Non-intersecting walks arose in yet another context, that of vertex 
models in statistical mechanics, where it was noticed that if the vertices of the six-vertex model 
are drawn in a particular way they could be interpreted as lattice paths |l(], [l(J. If one of the 
vertices had weight zero, giving a five-vertex model, the resulting paths were non-intersecting. 
The vertex models are traditionally solved by expressing the partition function (a generating 
function) in terms of transfer matrices. The partition function is then evaluated by either of two 
very powerful techniques, that of commuting transfer matrices Q or by direct diagonalisation 
of the transfer matrices using the coordinate Bethe Ansatz ||, O]. In this paper we bring 
the independent results of the two communities together for the case of N non-intersecting 
paths. We will show that the Bethe Ansatz (from statistical mechanics) and the Gessel- Viennot 
Theorem (from combinatorics) are essentially equivalent for a fairly general problem on the 
square lattice. The theorems proved should be easily generalisable to other planar lattices. The 
connection between the six-vertex model and non-intersecting path problems, from which this 



correspondence stems, has been recently discussed in general terms [10]. Here we shall consider a 
model equivalent to a generalised five-vertex model on the square lattice where the (Boltzmann) 
weights associated with walk edges are inhomogeneous in one direction. 



2 The model 

A lattice path or walk in this paper is a walk on a square lattice rotated 45° which has steps in 
only the north-east or south-east directions, and with sites labelled (m, y) (see figure 1). A set of 
walks is non-intersecting if they have no sites in common. We are concerned with enumerating 
the number of configurations of N non-intersecting walks, starting and ending at given positions, 
in various geometries: 1) walks in a plane without boundaries; 2) walks which are confined to 
the upper half plane; and 3) walks which are confined to a strip of a given width, L. More 
generally, one may be interested in interacting cases where the walks nearest the boundaries are 
attracted or repulsed by contact interactions: combinatorially this requires knowledge of the 
number of walks with particular numbers of contacts with each of the boundaries. In this paper 



we shall focus on case 3 since the other 2 cases can be easily obtained from this case as limits. 
To more easily describe our model we require the following sub-domains of Ij N 

Sl = {y 1 1 < V <L,y G Z and y odd}, (2.1a) 

Sl = {y | < y < L, y <E Z and y even}, (2.1b) 

SL = {y\0<y<L,yeZ}, (2.1c) 

Ul = {(yi,-.. , 2/jv) 1 1 <yi < ■■■ <VN<L, yi£S L } (2.1d) 

Ul = {(in,... ,2/jv) |0 < j/i < ... < 2/7v < L, yi^Sh] (2.1e) 

= {(yi,-.. ,1/jv) 10 < yi < ... < 2/7V < L, yi£S L } (2.1f) 

p o e 

We will use Ul to denote or Let A?~ non- intersecting walks, confined to a strip of 

p 

width L, start at y-coordinates y* = {y\, . . . ,y % N ) £ Ul in column m = of the lattice sites 
and terminate after t steps at y-coordinates = (y(, . . . ,ytr) € Ul in the column. If f 
is even then p' = p else p' = p, where p is the opposite parity to p. We will only consider 

e o / — IT I 1 \\ 

the case that L is odd so that | Ul \ = I Ul \ = ( 2 N )■ (If £ is even a null space enters 
the subsequent analysis of the transfer matrices leading to a distracting complication.) We are 
considering paths such that a) if (m — 1, y) is the position of a path in column m — 1 the only 
possible positions for that path in column m are (m, y') with y' = y ± 1 and < y' < L and b) 
the non-intersection is defined through the constraint that if there are ./V sites occupied at m = 
then in each column of sites (0 < m < t) there are exactly N occupied sites. We generalise the 
walk problem associated with the five- vertex problem [16, [l(]] by assigning a weight w(y,y') to 



the lattice edge from site (m — 1, y) to (m, y') with y' = y db 1 (see figure 1). Notice that, since 
w(y, y') is assumed independent of the column index m, due to the square lattice structure the 

P P 

weights are periodic in the m direction with period two: Note if y G5l then y €Sl, and in 
general w(y, y') ^ w(y', y). For the sake of generality we also associate an arbitrary weight v{y l ) 
with each of the sites occupied at m = 0. The weight associated with a given set of walks is the 
product of w weights over all edges occupied by the walks multiplied by the product of the v 
weights for each of the initial sites occupied. The generating function Zt (y* - *Y ), of N walks 
of length t starting at y = y* in column m = and finishing at y = in column m = t is the 
sum of these weights over all sets of walks connecting y l and y^ : 

=Af N 1 

Z t (y^yf) = J2 Y[v(Vj(0)) H^i™- VM™)) ( 2 - 2 ) 

y j=l m=l 

where yj(m) is the position of the j th walk in column m and the set y is given by 

y = {yj(m)\l <j<N,0<m<t,l<y 1 {m)< y 2 {m) < ■■■ < y N (m) < L, 

yj (m) = Vj (m - 1) ± 1 and Vj (0) = y), Vj (t) = y/.} (2.3) 

With homogeneous weights away from the boundaries but extra weights at the boundaries the 



associated six- vertex model has been considered in [051. However, we note that in [15|, and in 



most other studies of the six-vertex model, only such properties of the model are calculated that 
are averages over all numbers of walks, N. Here in contrast we are considering the generating 
function for a fixed number of walks, N, of a fixed finite length t. We will use the transfer matrix 
method from statistical mechanics to find this generating function in terms of a determinant of 
one- walk generating functions. 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 



Figure 1: Three non-intersecting directed walks of length t = 14 in a strip of width L = 
9. The variables m, y % - and y? shown. The walk closest to the lower wall has weight 



v(l)w(l, 2)w(2, 3)id(3, 2) • • • u>(l, 0)w(0, 1) with y\ = 1 and y{ 



f 



1. 



Precisely, we will show that the partition or generating function, Zt (y* — ► y^), for non- 
intersecting walk configurations of N walks in which the j th walk starts at j/j and arrives at 
after t steps is given by the following determinant: 



z t (y l 
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Z S t {y\ 
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(2.4) 



■ f 

where Zt {Vj — is th e generating function for configurations of a single walk starting at y*- 
and ending at y^ in a strip of width L. This is then a generalisation of the 'master formulae' 



of Fisher (equation (5.9) of Q) and of Forrester (equation (4) of ||) where unweighted non- 
intersecting walks are considered. Importantly, this determinantal result is precisely that given 
by the very general Gessel-Viennot Theorem & |8| for this problem: the walks considered on 
this lattice with the generalised weights considered satisfy the conditions of the theorem (see 



10] for a discussion of a one wall case with homogeneous weights). 



3 From Bethe Ansatz to determinant 

3.1 Transfer matrix formulation 

The generating function of our walk problem in a strip can be formulated as the matrix element 
of a product of so-called transfer matrices (see below for their precise definition). This 'transfer 
matrix' contains the weights of all the possible edge configurations of N walks of two adjacent 
columns of sites. As such it is related to an invariant subspace of the full six- vertex transfer 
matrix [pL0|J . In other words the non-zero elements of the six- vertex transfer matrix occur in 
diagonal blocks each of which is an "iV-walk transfer matrix" for some N. There are of course 
several ways of setting up a strip transfer matrix. The most common are the site-to-site and the 
edge-to-edge matrices which add either edges or sites respectively to the walks. In the case of 
the six-vertex model the edge-to-edge matrix is usually chosen. Our choice of vertex-to-vertex 
(we shall use vertex and site interchangeably) matrix is determined by the requirements that the 
corresponding eigenvalue problem be as simple as possible while ensuring the non-intersecting 
constraint can be easily implemented. 

The calculation of the generating function using the transfer matrix formulation then hinges 
on the spectral decomposition of the matrix. The matrices we diagonalise will be matrices that 
add two-steps to the evaluation of the generating function at a time. However the generating 
function is initially constructed in terms of one-step transfer matrices as these are useful since 
the non- intersecting condition is explicit. 

e o 

Definition 1. Let y G Sl and y' G Sl- For N = 1 the one-step transfer matrices are defined as 



CO 



, if \y - y'\ > 1 

Tl) , = { (3.5a) 

w(y,y') if|y-y'l = i 



y,y 



and 



,0 if W - y\ > 1 

Ti) , = (3.5b) 
™ (w(y',y) if|y-y'| = l. 

The iV-walk transfer matrices for N > 1 are constructed from sub-matrices of a direct product 



of the above N = 1 matrices: 

/ N \ 

oe I ^-^v oe \ o e 

(Tjv)y,y' = ( ® Ti I y € IAl an d y £ Ul (3.6a) 

Vi=i / y y 

and 

/ JV \ 

eo / eo \ e o 

(Tjv) y ',y = I ® Ti J y' eU L and y eU L (3.6b) 

Vi=l /y',y 

oe eo 

Note that since we are only considering L being odd we have that Tat and Ttv are square 



matrices and that in general ( Ti ) 7M Ti ) ■ More importantly, the restriction of the row 

\ Jy',y \ ) y,y> 

and column spaces of the direct product eliminates the possibility of two walks arriving at the 
same lattice point. Furthermore the condition that the one- walker transfer matrix vanishes for 
\y' — y\ > 1 prevents the generation of configurations in which pairs of walks "cross" without 
sharing a common lattice site (only nearest neighbour steps are allowed in all cases) . This "non- 
crossing" condition is unnecessarily restrictive in the one-walk case. However, for N > 1, if 
further neighbour steps are allowed then, in general, it is not possible to use the Bethe Ansatz. 

This condition is the analogue of the "non-crossing condition" of the Gessel-Viennot Theorem. 

=M . 

The generating function Zt (y l ~ * y ) of N non-intersecting walks of length t in a strip is 

=M 

related to Z 4 _i(y*— >y) by recurrence, the coefficients of which are the elements of one of the 
two one-step transfer matrices defined above. This relationship is given by the following lemma. 

Lemma 1. The generating function Zt (y* - * y ) is given for t > 0, depending on whether 
y J £Ul or not, by 



=M eo , o 

E _,<= -Zt-l(y l ^y)(Tjv) v ,y/ for yf eU L 



(3.7) 



oe 



E e ° Z t -i(y^y)(T N ) y , y f for y G Ul 



Proof. A simple proof of this Lemma can be constructed using induction on t. □ 
Together with the initial condition 

z"(y l ^y f ) = 5 yKyf V(y l ), (3.8) 



where 



N 



V(y i ) = l[v(yi), (3.9) 



a=l 



=N . 

equation ( p.7j ) determines Z t {y % — > y ) ■ 

A simple corollary of this Lemma (again shown by induction) is that the partition function 
can be written in terms of "two-step" transfer matrices. 



Corollary. The generating function Zt (y* — > y-* ) is given, depending on whether t is even or 
odd, as 



z 2r (y l ^y f ) 



y\y f 



for y l E U L and y f G Ul 



(3.10a) 



V{y l )[(T N ) r ) . , for y* £U L and yf €Ul 
' y\y f 



or 



=N . , 



(e-e \ eo 
(TK) '),., y 



WJXL:,; (T» 

y£Hi \ / y',y 



(Tjv)y, y / M y 4 G^l andy^GUL 

oe o „ e 

(Tiv) y , y i /or y 4 eZ^L andyt £Ul 



for r = 0, 1, 2, . . . , where the two-step transfer matrices are defined as 



e-e eo oe 

Ttv = TV Ttv 

o-o oe eo 

= Tn Tn- 



(3.10b) 

(3.11a) 
(3.11b) 



The two-step transfer matrices correspond to adding two steps to the paths. We show below 
that the generating function Zt {y l ~ > y ) may in fact be expressed in terms of the eigenvalues 
and eigenvectors of the two-step matrices. In general the two-step matrices are not symmetric 
and hence one has to consider left and right eigenvectors. The iV-walk eigenvectors will be 
constructed from the one- walk eigenvectors via the Bethe Ansatz. Used in this context the 
Ansatz expresses the components of the A^-walk eigenvectors as a determinant of the components 
of one- walk eigenvectors (see ( |3.17| ) below). 

3.2 From transfer matrices to determinants 



This section contains our main results in the form of two theorems that state under what 
conditions the A^-walk generating function can be written as the determinant fl2.4|)- In summary, 
our theorems proven below show that the equivalence of the Bethe Ansatz in the form of equation 
( 3,17| ) and the result of the Gessel-Viennot Theorem in the form ( p. 33 ) rests on showing that 
for any given problem the Ansatz is sufficiently good to provide a spanning (or complete) set 
of eigenvectors. This in turn depends only on the completeness of the one-walk eigenvectors 
(conditions of Lemma ||). 

In this paper it is not our purpose to prove that all choices of the functions w(y,y') which 
make up the elements of the one-walk transfer matrices allow for the conditions of our theorems 
to be satisfied. However, these theorems have been subsequently used to analyse some interesting 
cases in detail |J where as a consequence we show that an arbitrary boundary weight on one 
boundary (and special weight on the other) with homogeneous weights otherwise satisfies the 



conditions of the theorems. Also, we point out that since the Gessel-Viennot determinant holds 
for any function w(y, y 1 ) it is almost certainly true that our conditions are satisfied always. It is 
however our main purpose here to demonstrate how the Bethe Ansatz gives the solution of the 
iV-walk problem given the solution of the one-walk problem. 

The first theorem states the conditions under which the iV-walk transfer matrices can be 
diagonalised using a Bethe Ansatz: the major condition is that the one-walk transfer matrix 
problem can be solved — see Lemma pi The second theorem basically states that if the Bethe 
Ansatz gives a complete set of eigenvectors for the iV-walk problem for any N, then the iV-walk 
generating function is a determinant of one- walk generating functions. 

To begin we set up the conditions that define a solution of the one-walk problem: these will 
be the conditions our two theorems require. 

o 

Lemma 2. Suppose that there exist linearly independent sets of column vectors xfiufkeKi an d 

e 

{^k}keKu where K\ is some index set, which satisfy 

eo o e oe e o 

T\V% = \k<Pk and Ti<P% = \ k (3-12) 

eo oe 

with Afc € C, and which span the column spaces of Ti and Ti respectively (in which case they 

e-e o-o . . 

are said to be complete). Further let Ti and Ti be defined by (|3.1l| ) then 

o e o-o e-e 

(i) and are right eigenvectors of Ti and Ti respectively with eigenvalue Xf. 

o e 

(ii) corresponding sets {^Pj^kelCi and {'pfykeKi of row vectors may be found such that 

k*- h = 8k,k> and Yl k(v) k*(v') = <W (3.13) 

keKi 

for each p G {e, o}, where the * denotes complex conjugation. Note that the vectors of 



(3.13) have components indexed by y. 
(Hi) the row vectors of (ii) satisfy 

o oe e e eo o 

^Ti=A fe ^ and ^ Ti= A fe ^ (3.14) 

o e o-o e-e _ 

and also Vi and ^ are left eigenvectors of Ti and Ti respectively with eigenvalue Xf,. 

Note that since the vectors in ( 3.12[ ) span the space the cardinality of the index set /Ci is 
[L + l)/2. The proof of the lemma is elementary linear algebra and we omit it. Notice that if 

. . e e 

Afc is a solution of ( 3.12 ) then so is — with vector replaced by — Pf. These vectors are 
clearly not independent and normally sufficient independent vectors to form a spanning set are 
obtained by taking only the positive values of A^. 

From the above left and right one-walk vectors we now construct the A^-walk vectors and 

e-e o-o 

hence eigenvectors of T n and T n ■ 



e-e oo 

Theorem 1. Let Tat and Tn, N > 1, be given by equations ( |3.11| ) , By imposing an arbitrary 
ordering on the elements of K,\ define 



K-n = {k = (ki, &2, • • ■ k^)\ki € K,\ and k\ < k<z < ■ ■ ■ < k^} (3.15) 

and 



N 



a=l 



p 

(a) If for C € {L, R} and p G {e, o], {fk}k€Ki satisfy the conditions of Lemma [| then the 
v 

vectors {^}keic N given by the Bethe Ansatz, 



N N 

fey) = E e - II fe>*) = E e - II Kb**) y zul, (3.i7) 

crgP]v a=l o"6-Pjv a=l 

where Pn is the set of N\ permutations of {1, 2, . . . , N}, a = (o"i, 02, . . . , ctn) £ Pn and e a is 
the signature of the permutation a, satisfy 

oe e o eo o e 

TV = A k $k and Ttv $k = A k *k- (3.18) 

p p 

(b) Moreover the conclusions of parts (i), (ii) and (Hi) of Lemma\Q hold with <p? replaced by <j>£, 

eo oe eo oe 

Ti and Ti replaced by Tn and Tn, fci replaced by K-n , and replaced by A^. 
The proofs of part of this theorem and Theorem ^ require the following result. 
Proposition 1. For k € Kn and y £ £/l iei 

AT AT 

^k(y) = E ^II^^ and *k(y) = E e - 11^^)- ( 3 - 19 ) 

^4ko lei 

AT 



/(k) = [] /(*«) (3.20) 



then 



a=l 



N 



/(k)*k(y)* k (y') = E e -II E ffaWk-iVaWMj ) (3-21) 



anc 



A 



E *k(y)*k'(y)= E e -H E ^(i/«)^>«) ( 3 - 22 ) 



Proof. 



N 



E /(k)$ k (y)^ k (y') = E e - E E trX{f(K)^ kTa (y a )^ kaa (y' a 

kG/Cjv o-ePjv ke/CjvrGPiv Qf=l 



: E E E II f^^k^iy^kr^y'^ 

(t'&Pn ke/Cjv tSPjv a=l 



(3.23) 
(3.24) 



The double sum over permutations, t and k E /Cat is equivalent to summing each k a indepen- 
dently over /Ci (terms for which two or more components of k are equal make zero contribution) 
and the first result follows. The second result follows in the same way by interchanging the roles 
of k and y. □ 



Proof, (of Theorem |J) (a) We first obtain the cyclic property ( 3.18 ) as follows 



(tn Ik 



N 



e e e - nsu* 



y'eUL 



a=l 



N 



E ^ E II( T 0. v ( using O) 

o-ePjv 



E 

<tGPjv 



E (5.) 

,2/ieSi 



E (? 



A6 5i 



vn,v n 



o-ePiv 

A k ik(y) 



(3.25a) 
(3.25b) 

(3.25c) 

(3.25d) 
(3.25e) 



The critical step, and the whole reason for introducing the Bethe Ansatz, is to enable one to go 
from the restricted sums of (3.25b) to the unrestricted sums in ( |3.25c ). This is justified for two 
reasons, 

e e 

1. since 3?k(y') is a determinant, if any of the y a 's are equal then = - this allows the 
restriction y[ < y' 2 ■ ■ ■ < y' N on the sum to be relaxed to y[ < y' 2 ■ ■ ■ < y' N 

2. the y a are in strictly increasing order combined with the fact that the matrix elements of 

oe 

Ti, are only non-zero if \y a — y a \ < 1 allows the restriction on the sum to be removed 
altogether. 

The second part of ( |3. 18|) follows mutatis mutandis. 

e e-e 

(b) (i) The vector is a right eigenvector of Tat with eigenvalue A£, since 

e-e e eo oe e eo o „ e 

T7v$k = Tat Tat = T/vA k = K $k 



. U O'O 

which follows from ( |3.18| ). Similarly is a right eigenvector of Tat with eigenvalue A£. 



(ii) Let us start by deriving the first result of part (ii), namely the orthogonality and normali- 
sation condition. Using (3.22), for k, k' £ /Cat 



AT / \ 

E ^*(y) k>(y) = E e - II E feet**) 

p 



V 

yeWL 



o-ePjv a=l \ 



(using (gig )) 



= e ^ n s **M a 

ctGPjv a=l 
AT 

= n Sk «< k ' a 

a=l 

since the components of k and k' are in the same order only the identity permutation gives a 
non-zero delta function product. Thus 



p p 

$k • $k = °k,k' • 



(3.26) 



Our derivation of the second part of (ii), namely the "completeness condition", closely paral- 

p 

lels the derivation of the orthogonality and normalisation condition. Using (p. 21), for y, y GUl 



N 



E *S(y) #V) = E ^ II E KM. 



V 



£-114 

(TGPiV 

AT 

n* 



(using ( p3D ) 



AT 



so 



E fey)&V) = * 



y.y 



(3.27) 



keCjv 



Notice that \K N \ = ( 2 ^^) which is the row (and column) space dimension, as it should be for 
completeness. 

(iii) Using basic linear algebra gives 



and 



$kTjv =$k A k 



$k TV=$fcA£ 



e eo o 



o, oo o 



(3.28) 

(3.29) 
□ 



Lemma 3. If the conditions of Theorem [J hold th 



en 



zAy i ^y f ) = V(y i )J2 fey^L feV) y l eu L and yf eu L (3.30) 

ke/Cjv 

where if t is even, p' = p but otherwise p and p' are of opposite parity. 



Proof. If the conditions of Theorem @ hold then we have that (HH), (ET271) , &2§j ) and (^29|) 
are valid. Using (gj2§) and ( ggg ) it follows that 



(TjvW = E A k k(v) &V) and (Tiv) y ,y = E A k fe(y) &*(y')- (3-31) 



Also, using ( pT29| ) and flO^D it follows that 



(Tiv) y> y = E A k &(y) &V) and ( Tiv) y> y = E A * &(y) (3-32) 

kG^Cjv ke/Cjv 

Substituting these into ( |3. 1Q| ) and using ( 3.26| ) gives the result immediately. □ 



Theorem 2. If the conditions of Theorem [7] hold th 



en 



ZAy^y f ) = det \\Z S t {y i a -+V f f) )\\ 0h f>=i...N- (3-33) 



Proof. Using (3.9), ( 3.21| ) and ( |3.30| ) (which follows from ( p. 13 ) by Lemma 



-N N 



z t (y l ^y f )= E ^11 E *knM)rt*M 

ctGPjv a=l y k a eKi 
N =s 

= E ^\{z t {yi^yL)- (3-34) 

o-GPjv a=l 

which is an expansion of the required determinant. □ 

Finally, we point out that the result of Theorem || is also the conclusion of the Gessel-Viennot 
Theorem. 

3.3 One wall and no wall geometries 

When L > (t — \y l N — y^ N \)/2 + max(y]y, y^ N ) the walk closest to the wall at y = L cannot touch 

it since the N th walk (of t steps) needs at least \y l N — y^ N \ steps to go from y l N to y^ N and any 

excursion close to the wall from the end point closest to the wall (max.(y l N , y N )) needs just as 

many steps to return (hence the factor of 1/2). Hence, when this condition holds, the strip 
=s . , 

generating function, Z t {y % a ~^Up) becomes equal to the generating function for walks that are 
affected by only one wall, Zf(y l a — > y^). Hence taking the limit L —> oo gives the following 
corollary: 



P r P 

Corollary. For y l GUl and y> £Ul, the N-walk generating function with only one wall at 
height y = is given by, 

Zf(y i ^yf) = det\\Z?(yi^yf ) )\\ a , f3=1 ... N . (3.35) 

If we also condition the walk closest to the wall at y = so that it cannot touch that wall 
we will end up with the "no boundary" results. 
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